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We investigate two kinds of coreless vortices with axisymmetric and nonaxisymmetric configu- 
rations in rotating two-component Bose-Einstein condensates. Starting from the Gross-Pitaevskii 
energy functional in a rotating frame, we derive a nonlinear sigma model generalized to the two- 
component condensates. In terms of a pseudospin representation, an axisymmetric vortex and a 
nonaxisymmetric one correspond to spin textures referred to as a "skyrmion" and a "meron-pair" , 
respectively. A variational method is used to investigate the dependence of the sizes of the stable 
spin textures on system parameters, and the optimized variational function is found to reproduce 
well the numerical solution. In the SU(2) symmetric case, the optimal skyrmion and meron-pair are 
degenerate and transform to each other by a rotation of the pseudospin. An external rf-field that 
couples coherently the hyperfine states of two components breaks the degeneracy in favor of the 
meron-pair texture due to an effective transverse pseudomagnetic field. The difference between the 
intracomponent and intercomponent interactions yields a longitudinal pseudomagnetic field and a 
ferromagnetic or an antiferromagnetic pseudospin interaction, leading to a meron-pair texture with 
an anisotropic distribution of vorticity. 

PACS numbers: 03.75.Lm, 03.75.Mn, 05.30.Jp 



I. INTRODUCTION 



Since the experimental realization of quantized vor- 
tices in alkali atomic Bose-Einstein condensates (BECs) 
Qil^l^S 1 there have been a growing interest in new phe- 
nomena related to vortices in rotating BECs. A strongly 
correlated quantum Hall-like phase may appear in sys- 
tems which rotate so rapidly that the size of the vortex 
cores becomes comparable with the inter-vortex separa- 
tion "^1 . Another direction that has not yet been explored 
so much concerns rich vort ex p hases in multicomponent 

BECs iHiiiiaiiiiiiiiiiiiiiiiiiiiiiiaiiiii^ 

Since alkali atoms have hyperfine spin, multicomponent 
BECs can be realized if more than one hyperfine-spin 
state is populated in the same trap [2^ 0] . The systems 
described by the multicomponent order parameters al- 
low the excitation of exotic topological defects that have 
no analogue in systems with a single-component order 
parameter. For example, while a quantized vortex in a 
single-component order parameter should have a singular 
core, it is possible to excite a "coreless" (non-singular) 
vortex in multicomponent systems. On the analogy of 
such topological defects found in other physical systems 
such as superfluid ^He '2^, unconventional superconduc- 
tors quantum Hall systems jl^l, nonhnear optics 
psf. nuclear physics [l^ and cosmology js^l, creating 
exotic topological defects in atomic BECs could give us 
a further insight into related problems and offer a good 
opportunity to study their physical properties in much 
greater detail. 



This paper addresses the structure and the ener- 
getic stability of two kinds of coreless vortices in two- 
component BECs based on a generalized nonlinear sigma 
model (NLcrM). One of them has an axisymmetric struc- 
ture, in which the core of one circulating component is 
filled with the other nonrotating component. This vor- 
tex state was created by Matthews et al. l], where they 
utilized a phase imprinting technique by controlling inter- 
conversion between two components spatially and tem- 
porally with an external coupling field [31 . The structure 
and dynamical stability 10] of this vortex state have 
been studied. Here, the spinor nature of the order pa- 
rameter allows us to interpret the vortex state of the 
two-component BECs in terms of a "pseudospin" . The 
pseudospin representation of the two-component BECs 
reveals that the vortex state observed in Ref. cor- 
responds to a spin texture referred to as a "skyrmion" 
15, 29] or an "Anderson- Toulouse vortex" jOj, i31j. A 
spin-1 analogue of this texture was created by Leanhardt 
et al. ^] and theoretically studied in Ref. [l^, [iS ■ 

Another coreless vortex which we discuss in this paper 
has a nonaxisymmetric structure, which may be regarded 
as a pair of "merons" 0, |2^ |23] or "Mermin-Ho vor- 
tices" js^] • This configuration can be realized when each 
component has one off-centered vortex. Our previous 
study showed that an external driving field that couples 
coherently the internal hyperfine states of the two com- 
ponents stabilizes a nonaxisymmetric vortex state ^16j. 
Then, the internal coupling induces an attractive inter- 
action between the vortex in one component and that in 



2 



the other, forming a "vortex-antivortex molecule" , which 
is bound by a domain wall (branch cut) in the relative 
phase space . In contrast to a vortex-antivortex pair 
in conventional superfluid systems, this pair has circula- 
tions of the same sign in each individual phase space, but 
has the opposite sign in the relative phase space. When 
more than one vortex molecule is present, one component 
accomodates vortices and the other accomodates antivor- 
tices. 

It is known that the system of two-component BECs 
such as those studied by a JILA group P, approxi- 
mately possesses the SU(2) symmetry owing to the near- 
equal scattering lengths within and between the \F = 
l,mF — ~1) and \F = 2,mp = 1) hyperfine states of 
^"^Rb atoms. Here, the two kinds of coreless vortex states 
are degenerate in a completely SU(2)-symmetric system. 
However, adding an external perturbation that breaks 
the SU(2) symmetry would make the stability problem 
of those vortex states nontrivial. We investigate their 
structure and stability by exploring the NLcrM of the two- 
component BECs, which is derived from a pseudospin 
representation of the Gross-Pitaevskii energy functional. 
Using appropriate trial functions for the skyrmion or the 
meron-pair, we analytically obtain an almost exact spin 
profile for each case, which leads to a great improvement 
upon the Thomas- Fermi approximation 6]. While our 
previous paper used a solution of a classical NLctM 
directly as a variational function, a more general ansatz 
used here gives accurate profiles of the vortex states. 
The SU(2) symmetry is broken by the difference between 
the intracomponent and intercomponent two-body inter- 
actions, and by an external field which couples coher- 
ently the internal hyperfine states of the two components 
P, 13 ■ Under the pseudospin picture, they give rise to 
pseudomagnetic fields and a pseudospin-pseudospin in- 
teraction, either of which has a great influence on the 
structure of spin textures. We do not discuss compli- 
cated skyrmion excitations with topologically nontrivial 
spin profiles such as those studied in Ref . [13, llal • 

After formulating our problem in Sec. ^ we derive 
in Sec. IIIII the NLcrM that describes the two-component 
BECs. Based on the NLcrM, we determine the optimized 
structure of an axisymmetric skyrmion using a varia- 
tional method in Sec. IIVI We determine the optimized 
structure of a nonaxisymmetric meron-pair in Sec. In 
Sec. IVII we address the effect of axisymmetry-breaking 
contributions on the two kinds of spin textures. We con- 
clude this paper in Sec. IVIII 



thermore, their internal states are coupled coherently by 
an external driving field j34|. Viewed from the frame of 
reference corotating with the trap potential, the Gross- 
Pitaevskii energy functional of our problem reads 
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^17V^j|vI/,(r)p-f||vI/,(r)|^ 
+gi2|vI'i(r)nvE-2(r)p 
-TiCR [vI/*(r)vE'i(r)e-'^* + ^\{v)^^{v)e^'^'A , (1) 



where $7 = fiz, ^E*! and 5*2 denote the condensate wave 
functions in the two hyperfine states, and gi, (72, and 1712 
characterize the atom-atom interactions. Here gi, 172, 
and g\2 are expressed in terms of the s-wave scattering 
lengths a\ and a2 between atoms in the same hyperfine 
states and a\i between atoms in different hyperfine states 
as 



5i 



m 



(■ 1 47rn ai2 , , 
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The last two terms in Eq. describe a coherent cou- 
pling induced by an external driving field, which allows 
atoms to change their internal state coherently js^l • Since 
the driving field is time-dependent, we have introduced 
the frame in which the driving field is time-independent 
(i.e., the frame of a laser field). Here, a.'R(> 0) is the 
Rabi frequency and A is a detuning parameter between 
the external field and the atomic transition. Throughout 
this paper, we set A = for simplicity by assuming a 
complete resonance |35j . 

It is convenient to measure the length, time and en- 
ergy scale in units of 6ho = yjfi/mio., and Tilo, re- 
spectively. Renormalizing the wave function as ^'i 

%/iV*j/6^^^ with the total particle number N ^ N1+N2, 
and the energy as E/TlluN — > E, we obtain 
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II. FORMULATION OF THE PROBLEM 

We consider two-component BECs that are condensed 
into two different hyperfine states |I) and |2) such as 
those of *^Rb atoms. The two-component BECs are as- 
sumed to be trapped in the same harmonic potential 
V{y) = m{u!'^r'^ + Lo1z'^)/2. The potential is assumed to 
rotate at a rotation frequency about the z axis. Fur- 



Here, we denote the trapping potential as V = {(f — 
ri^)r^ 4- ci;^z^}/2 with a ~ uj^Ilo, the coupling constants 
as Ui = AnaiN/bho and U12 = 47rai2iV/6ho- Since the 
particles of one component can convert into the other 
if the internal coherent coupling is present, the total 
particle number N = Ni + N2 is conserved, and the 
normalization of the wave functions can be taken as 
/drd^'ip -I- |^'2p) = 1. Minimizing Eq. © with re- 
spect to ^'i and ^'2, we obtain time-independent coupled 



3 



Gross-Pitaevskii equations 
1 /V 



2 V i 



1 /V 
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-flxrj 'Si + y^fi + 

+ui2|«'2|^^'i -t^R*2 = M^'i, (4a) 

2 



-Ul2|*lP*2 -t^R*! M*2- (4b) 



Here, the chemical potential /i, which is common for both 
components, is determined by the normalization condi- 
tion. 



III. A GENERALIZED NONLINEAR SIGMA 
MODEL 



The pseudospin representation of the order param- 
eter with internal degrees of freedom is useful to ob- 
tain a physical interpretation by mapping the system 
to a magnetic system. Some physical properties of the 
double-layer quantum Hall system are well understood 
by projecting the system into a pseudospin space [2^ . 
Also, the spinor order parameter of two-component BECs 
allows us to analyze this system as a spin-1/2 BEG 
0,0,0,0,0,13. An exact mathematical correspon- 
dence can be established between these two systems. In 
this section, we derive the pseudospin representation of 
the energy functional, Eq. Q, which is the NLcrM that 
describes the two-component BECs. We assume that 
(^'2) corresponds to the up (down) component of the 
spin-1/2 spinor. The nonzero spin projection on the x-y 
plane implies a relative phase coherence between the up- 
and down-spin components. A similar NLctM was also 
discussed for a two-component Ginzburg-Landau energy 
functional l36|. 

We introduce a normalized complex- valued spinor x — 
[xi(r),X2(r)]^ = [|Xi|e'^S |X2|£^^]^ and decompose the 
wave function as '^i — -\/px(r)Xi(r), where pT is the 
total density and thus the spinor satisfies 



IXiP + IX2p 



(5) 



Substitution of the decomposed wave function into Eq. 
© yields 

E = y'dr|i(VV^)2 + ^[|VxiP + |Vx2p 

-2{n X r)(|xi|'V0i + |X2|'V02) + (O X r)^] 
+VpT - 2wrpt|xi||X2| cos(6'i - 62) 

+ \pt [co + ci(|xip-|X2|') + C2(|xip-|X2n']|, (6) 



where the new coupling constants are defined as 

U1+U2+ 2U12 
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(7a) 
(7b) 
(7c) 



The pseudospin density is defined as S = x(r)fx(r), 
where cr is the Pauli matrix. The explicit expressions of 
S = {S^,Sy,S^) are given by 

- (x^X2 + X2Xi) = 2|xi||X2|cos(0i -^2), (8a) 
^{xlX2 - Xki) = -2|xi||X2| sin(^i - 62), (8b) 
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5. = |xip-|X2p, (8c) 



where the modulus of the total spin is |S| = 1. 

Transformation of Eq. ((HJ to the pseudospin repre- 
sentation can be made as follows. On the kinetic-energy 
term, we have 
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and the second term of Eq. (jS)) becomes 
"(VS)2 
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Using Eq. Q and Eqs 
field 



(IXi|'V0i + |x2pV02-f^ xr) 

, we define an effective velocity 
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Sx^Sy), 



(9) 



which depends on the gradient of the total phase 8 
61 + 62 and that of the pseudospin. Thus, we obtain 



E= dv 



+VpT - ujrPtSx + ^(cq + CiS;, + C2SI) 

(10) 

This form is analogous to the classical NLcrM for Heisen- 
berg ferromagnets in which only the (VS)^ term appears 
[STf. In Eq. UlUI) . the four degrees of freedom of the 
original condensate wave functions and ^'2 (their am- 
plitudes and phases) are expressed in terms of the total 
density px, the total phase 0, and two of the spin den- 
sity components S = (Sx, Sy, Sz) (one of them is fixed 
because of |S| — 1) which are related to the relative den- 
sity and the relative phase. 

Unique features of Eq. (|1U|I that are absent in the clas- 
sical NLcrM, are: (i) the total density px, which is a pref- 
actor of the (VS)^ term and gives the pseudospin stiff- 
ness, is position-dependent because of the trapping po- 
tential, (ii) the total phase Q appears in the energy func- 
tional of Eq. (^UJ, (ui) the third term of Eq. 1(117}) gives 
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the hydrodynamic kinetic energy pT ( Voff — r2xr)^/2 asso- 
ciated with the topological excitation, and (iv) there are 
several anisotropic terms that break the SU(2) symmetry. 
The coherent coupling term with the Rabi frequency ljr 
works as a transverse (pseudo)magnctic field that aligns 
the spin along the x-axis. The interaction terms includ- 
ing the coefficients ci and C2 also break the SU(2) sym- 
metry. The coefficient ci can also be interpreted as a 
longitudinal (pseudo) magnetic field that aligns the spin 
along the z-axis. The term involving the coefficient ci 
determines the spin-spin interaction associated with Sz'^ 
it is antiferromagnctic for C2 > and ferromagnetic for 
C2 < [23. 

IV. AXISYMMETRIC SPIN TEXTURE: A 
SKYRMION 

Based on the NLcrM described by Eq. 1)10(1 . we dis- 
cuss the structure of vortex states in the two-component 
BECs. In the following sections, we consider the two- 
dimensional problem by assuming 3> w. Then, sepa- 
rating the degrees of freedom of the original wave func- 
tion as 5'i(r) = y)(j){z), we obtain the dimensionless 
two-dimensional GP equations 

/ y2 ^2 \ 
( — Y + y + "il"/^!!^ + ui2\H^ - flL.Ui = MiVdUa) 
/ y2 ^2 \ 

( ^ + y + U2 I V'2 P + Wl2 I V'l P - J ^2 = f^2^Uh) 

Here, we define effective two-dimensional coupling con- 
stants Ui = A'Kaii]N and M12 = A'Kai2riN with 77 = 
j dz\^{z)\V j dz\^{z)\\ 

A. Numerical results 

Figure n shows the density profile of the axisymmetric 
vortex state obtained by numerically solving Eq. Hll|l 
for ui = U2 = U12 = 1000 (cq — 1000 and ci = C2 = 
0), where the system possesses the SU(2) symmetry if 
= 0. Here, we assume that the ij^i component has 
one singly-quantized vortex at the center of the trap; we 
do not discuss the case in which the ipi component has 
a multiply quantized vortex. Then, the non-rotating 
component is located at the vortex core of the ■01 compo- 
nent; the core size is expanded due to the intercomponent 
repulsive interaction. As a result, the total density has 
no singularity and the condensates form a coreless vortex. 
This vortex structure was created by Matthews et al. Ij 
by following the phase imprinting method proposed by 
Williams and Holland (21 ■ 

The pseudospin texture corresponding to the axisym- 
metric vortex state is shown in Fig. I21 which was also 
discussed by Leonhardt and Volovik 9] and Mueller . 
At the center of the cloud, the "01 component vanishes, 
and the pseudospin points down in accordance with the 




FIG. 1: (Color online) (a) The density profile of the coreless 
vortex state consisting of the rotating xp\ component and the 
nonrotating i/)2 component for u\ = U2 = M12 = 1000 (co = 
1000, c\ = C2 uj-R = and Q, = 0.15. In this calculation, 

we fix the total particle number N = Ni + N2, but do not 
fix each particle number Ni. Then, the solution converges 
to iVi/iV2 = 2.465. (b) The cross sections of \ipif (sofid- 
curve), (dashed-curve) and the total density pT (dotted- 
curve) along the x-axis at y = 0, where bold and thin curves 
represent the results obtained from the numerical calculation 
and the variational calculation, respectively. 




FIG. 2; (Color online) (a) The pesudospin density S — X'^x 
corresponding to the vortex state in Fig. (b) The vectorial 
representation of the spin texture projected onto the x-y plane 
in the region [—6.7 < x,y < -1-6.7]. (c) The cross section of 
the spin texture along the x axis ai y — 0. 



definition of the spin Sz of Eq. The spin aligns 

with a hyperbolic distribution with (Sx,Sy) oc {x,—y) 
around the singularity at the center [Fig. |2l[b)]- At the 
edge of the cloud, the tp2 component vanishes, and the 
pseudospin points up. In between, the pseudospin rolls 
from down to up continuously as shown in Fig. |2Ic). 
This cross-disgyration spin texture is often referred to as 
a "skyrmion" in analogy to the work of Skyrme . 

Axisymmetric spin textures with continuous vorticity 
was extensively investigated in a field of superfluid '^He 
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B. Variational analysis 

To study the physical properties of the skyrmion in 
more detail, we make a variational analysis based on the 
NLcrM in Eq. The original NLcrM (Eq. ^ with 

only the (VS)^ term) admits a skyrmion solution and 
explicit analytic expressions are known (37l | . Here we 
take a more general form of the skyrmion solution as a 
variational function; the skyrmion solution of Fig. [21may 
be parametrized as |38j 



FIG. 3: (a) The distribution of the topological charge density 
g(r) of Eq. 11411 and (b) the vectorial plot of the effective 
velocity field Vcff of Eq. JSJ, corresponding to the solution 
shown in Fig. Q The insets represent the cross sections of g(r) 
and |vcff j along the y = line within the range —5 < a; < 5. 
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(15) 



|25| . When the condensate wave function is parametrized 
as 



5. = 



V'2 



(12) 



the configuration satisfying the boundary condition 
/3(0) — IT and /3(oo) = is referred to as an "Anderson- 
Toulouse (AT)" vortex [HsJl. We also have a "Mermin- 
Ho (MH)" vortex or a "meron" texture with the condition 
/3(0) = TT and /3{oo) = 7r/2 9, 32]. In the case of super- 
fluid '^He, a MH vortex is stabilized by the boundary 
condition imposed by a cylindrical vessel. However, in 
atomic-BEC system there is no constraint at the bound- 
ary; the value /3(r) at the boundary should be determined 
self-consistently as discussed later. 

It is known that the skyrmion has a topological invari- 
ant defined in a two-dimensional system as 
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(13) 



which is called a topological charge or the Pontryagian 
index The skyrmion with any spin profile is shown 
to have Q = ±1, whose sign depends on the direction of 
the circulation of a vortex. The integrand of Eq. H13|l is 
the topological charge density associated with the vortic- 
ity derived from the effective velocity Ves 

q{r) = -^e'-'S • d,S X - -^(V X Veff), (14) 

OTT ZTT 

where we used the relation J^i Xi'^Xi ~ ~ J2i Xi'^Xi 
{i — 1,2) in obtaining the last equality. The topological 
charge density q{r) characterizes the spatial distribution 
of the skyrmion. Figure 01 shows the spatial distribution 
of g(r) and the corresponding Veff-field. The topological 
charge is distributed around the center and, contrary to 
a conventional vortex in a single-component condensate, 
jvoff I vanishes at the center. This makes a coreless vortex 
without a density dip in the total density. 
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with |S| = 1. Equations H15|l with a = correspond to 
the explicit skyrmion solution of the classical NLcrM [s^l ■ 
The variational parameters A and a determine the size 
and the shape of the skyrmion. Typically, A represents 
the size of the domain in which the spin is reversed and a 
the degree of spatial variation of the spin inversion. If we 
take both A and a as variational parameters, the number 
of particles in each pseudospin component is not con- 
served (we will discuss below the situation in which each 
particle number is fixed). Therefore, the energetically 
stable configuration of the skyrmion brings about an im- 
balance in the particle number of each component. Such 
population imbalance can be controlled by the turning- 
off time of the coupling drive in the phase imprinting 
method 0, For the limit A — > (cxd) with fixed a 
all spins point up (down) along the z-axis, which means 
perfect polarization of the particle number as iVi 
{N2 N). Alternately, for fixed A, the limit A^i ^ A^ 
{N2 ~» A^) corresponds to a +00 (—00). Because 
the V'l component has one singly-quantized vortex at the 
center, the total phase Q in Eq. (|10|l is given by 



e = tan-^-. (16) 

X 

Substituting Eqs. (UHl and into Eq. yields 



E = dr< -(VVpt) + V'cffPT - WRPT' 
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(17) 



Now the degrees of freedom in the energy functional have 
been reduced to the total density and the two varia- 
tional parameters A and a. Here, we have introduced an 
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FIG. 4: (a) The eflFective potential Kft of Eq. ([THJ with a = 
(for r2 = and ajR = 0) for several values of A. A pure har- 
monic potential (obtained for A — > oo) is shown by a dotted 
curve. The inset shows the corresponding total density pro- 
files pt for Co = 1000. 
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FIG. 5: The total energy as a function of the skyrmion size 
A for several values of the rotation frequency Q. (fi =0, 0.05, 
0.1, 0.15) for Co = 10^ and ci = C2 = 0. The energy of 
the nonvortex state is 11.894. The inset shows the optimized 
value of a as a function of A. 



effective confining potential 

r2-^4A2e-"'-'[(ar2-hl)2-hl] 
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(18) 



which is radiaUy symmetric and determines the shape of 
the total density pT- Note that the term with the Rabi 
frequency wr, in Eq. (|17|l breaks the axisymmetry of the 
problem; we will discuss the effect of this term in Sec. Ivl 
and IVII Ignoring this term, we can calculate the total 
density pT by solving the equation 
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= /i, (19) 



where the chemical potential p is fixed by the normaliza- 
tion condition / drpx = 1. 

Before proceeding the analysis further, it is instruc- 
tive to investigate the form of the effective potential Vcs- 
For simplicity, let us first consider the case with a = 0, 
i.e., the classical skyrmion solution. Figure 0] shows the 
effective potential and the total density profile pt 
calculated numerically from Eq. (|19(l for several values 
of A for Co = 1000, ci = C2 = and = 0. For A = 
[N = Ni), the value of Vcs at r = diverges, which 
implies that the total density px should vanish at the 
center. This corresponds to the singular vortex core of 
the tpi component. With increasing A the central peak 
of VcS decreases and the total density becomes smooth 
because of the appearance of the nonrotating ip2 com- 
ponent, which fills the vortex core of the ipi component. 
The limit A ^ oo (iV N2) describes the vanishing i/^i 



component with vorticity. Then, the effective potential 
becomes a pure harmonic potential and the total density 
becomes |V'2p exactly. This character of Vcs also holds 
for a ^ 0. To describe px except for A ~ 0, therefore, 
we can employ the Thomas-Fermi approximation which 
neglects the quantum-pressure term — (V2.^/p^)/2y/px in 
Eq. (|17|l . Then, the total density is directly given from 

Eq. dm. 



1. SU (2) -symmetric condensates 

Using the condition J drpx = 1 , we numerically calcu- 
late the chemical potential in Eq. H19|) with the Thomas- 
Fermi approximation (the region with negative px is ne- 
glected) and the total energy i? as a function of A and 
a from the resulting total density px. Figure |21 shows 
the total energy E of the axisymmetric vortex state as 
a function of A for co = 10'^, ci = C2 = and several 
values of the rotation frequency Q, where the value of 
a is optimized for each A. For values of higher than 
0.17, more vortices enter the system. In the case with- 
out rotation (fi = 0), the energy decreases with A and 
suddenly drops to the value at which the system has no 
vorticity, i.e., complete polarization of the particle num- 
ber as N — N2 (see the inset of Fig. |5j , where there is 
no energy minimum corresponding to the vortex state for 
fixed A because a — > — 00. This means that the global 
minimum for = is a nonvortex state for any values of 
Ni and N2 under the fixed total particle number N, be- 
cause the interaction-energy terms now satisfy the SU(2) 
symmetry so that the energy of the nonvortex state is de- 
generate under the change of the ratio iVi /N2 with fixed 
N. However, beyond a certain critical value of fl, there 
appears an energy minimum as shown in Fig. [S] Since 
the minimized energy is lower than that of the nonvortex 
state at A — > -l-oo, the external rotation can make the 
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FIG. 6: The optimized value of the variational parameter A 
(solid curve) and that of a (dashed curve) as a function of 
rotation frequency Q for co — 10^, 2.5 x 10^, 10* and ci = 
C2 = 0. 
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FIG. 7: The bending angle /3(r) = cos~^ at the Thomas- 
Fermi boundary r = Rtf at which the total density pT van- 
ishes, as a function of Q. The optimized values of the varia- 
tional parameters A and a in Fig. [Sjare used. 



corresponding skyrmion texture stable globally. 

Figure|n|shows the values of A and a that give the min- 
imum of the total energy for cq = 10'^, 2.5 x 10^, and 10^. 
The size of the skyrmion decreases with increasing fi, as 
revealed by a decrease (an increase) in X„iin (amin)- The 
energy minimum appears at a certain critical frequency of 
n, where Amin (amin) diverges to +oo (— oo) and the crit- 
ical frequency decreases as cq increases. Therefore, the 
condensates with larger co can have a stable skyrmion 
at a lower rotation frequency. However, the full numeri- 
cal calculations of Eq. Hll|l show that additional vortices 
are nucleated and form a lattice of skyrmions beyond 
n ~ 0.17, 0.11, 0.05 for cq = 10^, 2.5 x 10^, 10*, re- 
spectively. These structures cannot be described by Eq. 
((T^ and the investigation of such periodic solutions re- 
mains to be made. From the comparison with the nu- 
merical solution as shown in Fig. H^b), the optimized 
variational functions Eqs. p5|l and 1)16(1 reproduce an 
almost exact numerical solution. Hence, our approach 
improves greatly the analytic treatment for the vortex 
states based on the usual Thomas-Fermi approximation 
studied in Ref. 1^. 

The bending angle (3{r) — cos^^ Sz in Eq. ifT^ de- 
creases smoothly from tt at r = as r increases. Figure 
shows the values of /3(r) evaluated at the Thomas-Fermi 
boundary r — i?TF at which the total density pT van- 
ishes, as a function of rotation frequency fl. The value 
/3(i?TF) can change arbitrary with f2 by varying the ra- 
tio N1/N2, which implies that an intermediate configu- 
ration between a MH vortex (/3(i?TF) = 7''/2) and an AT 
vortex (/3(i?TF) = 0) can be made thermodynamically 
stable. It is the ratio N1/N2 that determines (3{Rtf), 
which does not change for fixed N1/N2 as explained in 
Sec. IIV B 31 We also note that there exists a minimum 
value of P{Rtf) 0.09), below which additional vor- 
tices are nucleated. Therefore, an AT vortex can never be 
the ground state of the SU(2)-symmetric two-component 
condensates. A similar discussion is made in Ref. |l3 . [T^ 
for the case of ferromagnetic spin-1 BECs. 



2. ci and C2 dependence 

Next, we discuss the effects of ci and C2 in Eq. H17|l . 
the values of which are controllable by changing the scat- 
tering lengths. The coefhcient ci, proportional to the 
difference of the intracomponent interactions ui and 1*2, 
may be regarded as a (pseudo)magnetic field that aligns 
the spin along the z-axis; the positive (negative) ci aligns 
the spins downward (upward). Figure |S| shows the opti- 
mized value of A and a as a function of ci for f2 = 0.15 
and Co = 10'^. For the negative Ci side the size of a sta- 
ble skyrmion shrinks as \ci \ increases, which implies the 
spin alignment to upward with increasing a fraction of 
the rotating ■01 component. In this case, an AT vortex 
can be thermodynamically stable. Below A at ci ~ —85, 
the Thomas-Fermi approximation breaks down for small 
A owing to the appearance of the singular vortex core 
in pt- For the positive ci side the value of A increases 
rapidly and eventually goes to infinity (concurrently, a 
goes to —00), corresponding to the complete "spin-down" 
alignment (occupation of the nonrotating ip2 component 
only). This dependence on the sign of ci reflects the fact 
that the occupation of the component with vorticity costs 
a larger energy than that of the nonrotating component. 

The sign of C2 determines the nature of the spin-spin 
interaction associated with S^. For positive C2 (ui + U2 > 
2^12) the interaction is antiferromagnetic, preferring to 
the spatial mixture of the spin-up and spin-down compo- 
nents. While for negative C2 {ui+U2 < 2^12), the system 
enters the ferromagnetic phase, where the spin domains 
are spontaneously formed. This feature also appears in 
the stable structure of the skyrmion. Figure IHl shows 
the optimized values of A and a as a function of C2 for 

— 0.15, Co — 10'^ and ci — 0. For the antiferromagnetic 
case C2 > 0, we see no significant change of the spin pro- 
file compared with that of C2 = (the SU(2) symmetric 
case). For the ferromagnetic case C2 < 0, the variations 
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FIG. 8: (Color online) The optimized values of the variational 
parameters A (solid curve) and a (dashed curve) as a function 
of the longitudinal (pseudo)magnetic field ci for cq — 10^, 
C2 = 0, and S7 = 0.15. Insets show the cross sections of the 
variational function \ipi\'^ (solid curve), |V'2|'^ (dashed curve), 
and the total density pT (dotted curve) along the y = Q line 
(left: ci = -50, right: ci = 8). 
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FIG. 9: (Color online) The optimized values of the variational 
parameters A (solid curve) and a (dashed curve) as a function 
of the spin-spin interaction strength C2 for co ~ 10^, ci — 0, 
and f2 = 0.15. The insets show the cross sections of the 
variational function (solid curve), IV'al^ (dashed curve), 

and the total density px (dotted curve) along the y = Q line 
(left: C2 = -50, right: C2 = 50). 



of A and a are similar to that of FiglSJ where the size of 
the skyrmion shrinks with | C2 1 . Since each particle num- 
ber is not conserved in this calculation, the energetically 
favorable configuration in a ferrmagnetic phase tends to 
an overall spin polarization, i.e., complete polarization of 
the particle number. In this regime, there are two energy 
minima, one of which leads to the perfect polarization of 
the ^1 vortex state (corresponding to the results shown 
in Fig. O and the other leads to the ■02 nonvortex state 
(characterized by A — > oo and a — s- — oo; not shown), re- 
spectively, and which state possesses the global stability 
depends on the rotation frequency. 



Energy minimization under the fixed particle number of 
each component 



So far, we have allowed a change in the particle 
number of each component to calculate the minimized 
skyrmion size. However, the experiments at JILA on two- 
component BECs |23j were made under the condition in 
which each particle number is fixed. This restriction can 
be taken into account by noticing the relation 



N2 jdvpT{l-S,y 



(20) 



for a given A the value of a is uniquely determined by Eq. 
(|20|l . where both px and Sz are functions of two varia- 
tional parameters. Therefore, the energy minimization 
can be done with respect to one variational parameter 
which we choose to be A. As discussed before, this ra- 
tio determines the boundary value of the bending angle 
/3(r). 

Here, we consider the situation of the equal particle 
number iVi /N2 — 1 and investigate the stable size of the 
skyrmion as done before. In this case, we find that the 
stable size is not affected by the change of Vl. Figure [TUI 
shows the optimized values of A and a as a function of ci 
and C2. In this case, no complete polarization of the par- 
ticle number occur, so that the skyrmion may exist for 
all values of ci and C2; in this case, however, even if the 
optimal variational parameters exist, they do not ensure 
a local minimum of the total energy. In Fig. llOf a). for 
positive Ci the divergence of A„iin seen in Fig. |S1 is sup- 
pressed due to the conservation of each particle number. 
Alternately, to enlarge the "spin-down" domain where 
the "02 component is occupied, the size of the vortex core 
of the V'l component becomes large by increasing both 
Amin and aniin- For negative ci, the domain with the ro- 
tating 01 component tends to increase. As a result, the 
vortex core of the -01 component shrinks and a part of 
the "02 component is pushed out of the -01 component. 
Here, the total energy of the skyrmion with the same 
|ci| is always lower for positive ci than for negative ci. 
This is because it is favorable for the 0i component with 
the relation ui > U2 to spread spatially by the centrifu- 
gal effect with a vortex, and increase the "spin-down" 
domain. This is consistent with the experimental ob- 
servation that the vortex in the component having a 
larger intracomponent interaction is stable; otherwise it 
is unstable (see also the argument of the dynamical sta- 
bility in Ref. IQJ). In the C2 dependence, while we see no 
drastic change in the antiferromagnetic C2 > range, for 
C2 < a rapid increase of Amin and amin with increasing 
|c2| is seen in Fig. llOf b). This means that the spin-up 
or spin-down domains grow and their boundary becomes 
sharp with increasing |c2|. 
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FIG. 11: (Color online) (a) The density profile of the conden- 
sates Iv^il^ and IV'al^ for ui — U2 = ui2 — 1000 (co = 10'', 
ci = C2 = 0), f2 = 0.15 and ui-r = 0.05. The two components 
have the same particle number, (b) The cross sections of 
It/iiP (solid-curve), \tp2\^ (dashed-curve) and the total density 
Pt (dotted-curve) along the a:-axis at j/ = 0, where bold and 
thin curves represent the results obtained from the numerical 
calculation and those obtained from the variational calcula- 
tion, respectively, (c) The gray scale plot of the relative phase 



FIG. 10: (Color online) The optimized values of the varia- 
tional parameters A (solid curve) and a (dashed curve) as a 
function of (a) ci and (b) C2 under the condition of a fixed 
particle number of each component N1/N2 ~ 1 for co = 10'^, 
Q, = 0.15 (c2 = and ci = for (a) and (b), respectively). 
The insets show the cross sections of the variational functions 
IV^ip (solid curve), I'i/'al^ (dashed curve), and the total density 



Pt (dotted curve) along the y = line ((a) left: ci 
right: ci = 50 and (b) left: C2 = —30, right: C2 — 30) 



= -50, 



V. NONAXISYMMETRIC SPIN TEXTURE: A 
MERON-PAIR 

A. Numerical results 

In this section, we discuss a nonaxisymmetric vortex 
state. This structure can be regarded as a spin texture 
consisting of a pair of merons p7| or MH vortices [2Ml3^ . 
Our previous study showed that this meron-pair can be 
stabihzed thermodynamically in the presence of rotation 
and internal coherent couphng . A typical solution of 
this structure is shown in Fig. ^] Each component has 
one off-axis vortex and the density peak of one compo- 
nent is located at the vortex core of the other component. 
This results in a coreless vortex in which the total density 
has no singularity. 

Because of the presence of the coherent coupling, the 
profile of the relative phase ^(r) =61—62 plays an im- 
portant role in optimizing the structure, which is shown 
in Fig. Illf c). The relative phase shows that the central 
region is characterized by the configuration of a vortex- 



antivortex pair. In other words, the two vortices are con- 
nected by a domain wall of the relative phase with the 27r 
difference, which is characteristic of a topological soliton 
in two-component BECs with the internal coherent cou- 
pling. Son and Stephanov [s^l obtained an exact form 
of the domain wall in a homogeneous one-dimensional 
system. In our case, the one-dimensional profile of 0(r) 
along the x = line approximately takes a form 



my) 



-4 tan" 



-1 ^ky 



(21) 



where k^^ = (jV'i | |'(/'2|/2wr/Ct)"'^^^ gives a characteristic 
size of the domain wall, and the constant C = 2tt (0) for 
y > {< 0) makes the branch cut at y = shown in Fig. 
Illf c). Then, the vortex in one component and that in 
the other can be bound by this domain wall, forming a 
"vortex molecule" 

The pseudospin profiles corresponding to the solution 
of Fig. ^2 are shown in Fig. I12f a). While the profile 
of Sx is axisymmetric, Sy and Sz show dipole structures. 
The structure of S has a form similar to what is obtained 
by exchanging Sx and Sz in Fig. El The corresponding 
spin texture is shown in Fig. I12r b) and (c). The spins 
are oriented in the x direction everywhere except in the 
central domain-wall region, where they tumble rapidly 
by 27r. There exist two points corresponding to the lo- 
cations of vortices at which S is parallel to the z-axis. 
The spin around the singularity with S = +z (S = — z) 
has a radial (hyperbolic) distribution, characterized by 
{Sx,Sy) (X {—x,—y) (oc {x,—y)), and rotates through 
90° as it goes outward to become perpendicular to the x- 
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FIG. 12: (Color online) (a) The pesudospin density S = X'^X 
corresponding to the vortex state in Fig. 1111 (b) The vectorial 
representation of the spin texture projected onto the x-y plane 
in the region [—6.7 < x,y < 6.7]. (c) The cross section of the 
spin texture (b) along the x axis at y = 0. 



axis. This texture is known as a "radial- hyperbolic" pair 
of merons, which has been discussed in the study of topo- 
logical defects in superfluid '^He and a double-layer 
quantum Hall system ■ Since the trapping potential 
is axisymmetric, the energy is degenerate with respect 
to the orientation of the molecular polarization. If the 
molecule is polarized along thew-axis, the texture forms 
a "circular-hyperbolic" pair [I3I and can therefore con- 
tinuously change into the radial-hyperbolic form. 

We find that, when ui = U2 = U12 (ci ~ C2 ~ 0), the 
topological charge density q{r) and the effective velocity 
field Voff exhibit a radially isotropic profile as in Fig. 
even though each component forms a nonaxisymmetric 
vortex configuration. As we discuss later, this is due to 
the fact that the meron-pair texture and the axisymmet- 
ric skyrmion texture are equivalent as a topological exci- 
tation in the case of mi = U2 = U12] they can transform to 
each other through an overall rotation of the pseudospin. 
Therefore, the Rabi coupling alone does not break the 
axisymmetry of the topological excitation. However, in- 
clusion of both the Rabi coupling and unequal coupling 
constants ui 7^ U2 7^ U12 induces an anisotropy of the 
meron-pair as discussed in Sec. IVII 



B. Variational analysis 

In the case of ui = U2 = M12 (ci = C2 = 0), the spin 
profile of the meron-pair in Fig. 1121 is radially isotropic 



and may be parametrized as 



Sy 



(22) 



Here we assume that the vortex molecule is polarized 
along the x-axis. Compared with the ansatz solution in 
Eq. IjlSfl for a skyrmion, only Sx and Sz are exchanged. 
In this case, the ratio of the particle number given in Eq. 

is always unity, because Sz is an odd function. In 
Ref. we set a = to simplify the variatioal problem, 
and obtained a reasonable agreement with the numerical 
solution. The full variational analysis with two varia- 
tional parameters A and a gives almost perfect quantita- 
tive agreement with the numeircal one as seen below. 

The locations of the vortex cores are determined by 
two extremes of Sz, given by 



4A^ 



(23) 



It solution is represented as a; = ±M^(4aA^)/ y/a with the 
product log function W{z) [3^23l- Then, it is natural 
to take the form of the total phase 8 as 



= tan 



y 



X - 2Ae-"'~V2 



tan 



y 



(24) 



Substituting Eq. ^ and Eq. ^ into Eq. (jTUl), we 
obtain the total energy similar to Eq. H17|> : 



E 



dr 



Pt 



4A2 



r2 + AX^e^ 



(25) 



Here, the effective confining potential T4ff is the same as 
in Eq. (|18|l . Contrary to the axisymmetric skyrmion, the 
internal coupling term with does not break the radial 
symmetry of the problem. 

Because Eq. (|25f) gives the same energy for any value 
of A and a as that of the axisymmetric skyrmion given by 
Eq. I|17|l . if wr, = 0, i.e., for the SU(2) symmetric case, 
the skyrmion and the meron-pair have the same optimal 
values of A and a and their energies are degenerate. In 
other words, the energy is degenerate with respect to the 
rotation of the overall pseudospin between the two spin 
textures. We can see that turning on the Rabi term cjr 
always decreases the free energy of the meron-pair from 
any solutions of vortex states with the topological charge 
Q = +1 for cjR = 0. Therefore, adding the infinitesimal 
value of cijR is enough to stabilize the nonaxisymmetric 
configuration if ci = C2 = 0. 

Since the Rabi frequency may be regarded as a trans- 
verse magnetic field along the x-axis, the stable size of 
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i.e., the vortex molecule, was discussed previously |l6j |. 
For well-separated merons, the repulsive interaction be- 
tween them originates from the second and third terms of 
Eq. (|10|l . which are the gradient energy of the pseudospin 
and the hydrodynamic kinetic energy of the Voff-field. On 
the other hand, the attractive force between two merons 
are estimated from a tension of the domain wall of 
the relative phase to be T^d^; for a homogeneous sys- 
tem T(j = 8\tjji\'^\tjj2\'^k/pT with the characteristic domain 
size fc-i = (|^/'i||V'2|/2wrpt)'/' |33|, we have Td (x 
Then, the competition between the repulsive force and 
the attractive force creates an energy minimum so that 
the two vortices can form a bound pair. The contribu- 
tion of the other terms are almost constant except for 
small dm (< 0.30 for cq = 10'^), where the vortex core 
appears in the total density. Then, the Thomas-Fermi 
approximation cannot apply to the evaluation of the to- 
tal energy; Amin drops suddenly to zero with increasing 
because of the lacking of the energy barrier associated 
with the quantum pressure of px- The numerical result 
shows that the separation d^ decreases smoothly to zero 
with wr. 



FIG. 13: (Color online) (a) The optimized values of the vari- 
ational parameters A (solid curve) and a (dashed curve) as a 
function of lor for co = 10'', Q, — 0.15, ci — and C2 = 0. 
(b) The corresponding size of the meron-pair 2dm calculated 
from Eq. ^ (see Fig. ^h) for definition) for co = 10^ We 
also show 2dm obtained from numerical simulations as filled 
circles. The insets show the cross sections of the variational 
functions (solid curve), \ip2\'^ (dashed curve), and the 

total density pT (dotted curve) along the y = line (left: 
u-R. = 0.1, right: lor, — 0.5). 



the meron-pair as a function of ton. shows a similar be- 
havior as that shown in Fig. |H| The difference is that 
the Rabi term is proportional to instead of p^, which 
changes the behavior from that of Fig. |S1 In the presence 
of slow rotation '41*1, we calculate the minimized values 
of A and a as a function of wr, the result being shown 
in Fig. I13f a). As ojr increases, the minimized value of 
A (a) becomes smaller (larger), and eventually vanishes 
(diverges). This behavior corresponds to a decrease in 
the size dm of the meron-pair (see in Fig. Illf b) for the 
definition) as shown in Fig. Il^t b). where dm is calcu- 
lated from Eq. (|23|) . This indicates that the binding 
of the meron-pair becomes stronger with increasing wr , 
which is caused by the tension of a domain wall in the rel- 
ative phase between the two-component wave functions 
p^, "SjI . The variational result agrees perfectly with the 
numerical result. Beyond wr ~ 3.0 the separation dm 
vanishes, where the locations of the density nodes over- 
lap despite the intercomponent repulsive interaction. 

The physical origin of the binding of the meron-pair, 



C. Basis transformation of the pseudospin 



In the SU(2) symmetric case (ci = C2 = 0), a fur- 
ther insight of the nonaxisymmetric vortex is obtained 
by rotating the basis of the spinor so that the internal 
coupling becomes simpler. Note that the spin profile of a 
skyrmion [Fig. |2] and a meron -pair [Fig. 112) is connected 
only through the exchange of Sx and Sz ■ In this case, the 
Rabi term makes the cc-axis as a preferred one. Actually, 
if we rotate the spinors as V'± = (V'l i '^2)/V2 (the ba- 
sis along with the x-axis), the coupled Gross-Pitaevskii 
equations become 



1 /V 



2 V i 



$7 X r +V + cqPt 



V'+ = (m + ^r)V'-h, 



Here, the internal coupling is just the chemical poten- 
tial difference between the and "— " components. 
Then, the nonaxisymmetric structure in Fig. Illl is trans- 
formed to the axisymmetric vortex state of Fig. ^ i.e., a 
skyrmion. Here, the vortex core of the component 
is filled with the nonrotating "— " component. As one in- 
creases wr the number of the "— " particle drops because 
of a decrease in the chemical potential of tp- , and the vor- 
tex cores eventually become empty, which corresponds to 
the overlap of the vortex cores. 
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VI. EFFECTS OF SYMMETRY-BREAKING 
TERMS 



In this section, we discuss the effect of the symme- 
try breaking terms, which are neglected in the last two 
sections, by numerically solving Eqs. Hll|) . For an ax- 
isymmetric skyrmion texture, the Rabi term breaks the 
axisymmetry of the spin texture, and for a meron-pair 
the ci and C2 term cause symmetry breaking. 

First, we neglect the spin-spin interaction term (c2- 
term) and investigate the dependence of the stable struc- 
ture on ci and wr, which corresponds to a situation in 
which a longitudinal and transverse (pseudo)magnetic 
field are applied simultaneously. We first prepare a sta- 
ble axisymmetric skyrmion by applying the longitudinal 
magnetic field ci = —10, and then turn on the transverse 
magnetic field ui^. Figures [UT aV fdl show the variation 
of the spin texture as well as the cross section of the 
condensate densities along the y = line when the Rabi 
frequecy wr is increased. The transverse magnetic field 
shifts the skyrmion to the off-center and convert it into a 
meron. Correspondingly, another meron enters from out- 
ward and forms an asymmetric vortex molecule shown in 
Figs. I14f c') and (d). This change is a second-order transi- 
tion; there is no energy barrier to destabilize the axisym- 
metric skyrmion with respect to the transverse magnetic 
field. Figure IT^ el shows the topological charge density 
q{r) along the x-axis (the y = line). Although its dis- 
tribution is shifted from the center with increasing wr, 
there is no dramatic change in the global shape from the 
isotropic one. As wr increases further, the peak of q{r) 
gets back to the center and its value goes to infinity when 
the two merons merge. 

Next, we study the effects of the spin-spin interaction 
term (c2 term) on the stable structure of a meron-pair un- 
der a transverse magnetic field with wr = 0.05 (ci is fixed 
to be zero). Then, C2 changes the structure of the meron- 
pair dramatically, in contrast to the case of ci. Figure lT^ 
shows the equilibrium structure of the condensate den- 
sity, the spin texture and the topological charge density 
g(r) for the antiferromagnetic case (c2 = —20) and the 
ferromagnetic case (c2 = 20). For the antiferromagnetic 
case there is no significant difference in the density and 
spin profile, compared with the solution of C2 = in Figs. 
II II and [T^ However, the topological charge density, i.e., 
vorticity, is distributed anisotropically in such a manner 
that the distribution is elongated along the direction of 
polarization of the meron-pair. For the ferromagnetic 
case, the spin domains are formed, which gives rise to 
a considerable change in the density and spin profile as 
seen in Fig. llSf b). Most of the spins align up or down 
along the z-axis on both side, sandwiching the domain 
wall across which the spin flips rapidly. If the Rabi fre- 
quency o^R is increased further, the spins on both sides 
are laid along the a;-axis, converting into the well-defined 
meron-pair as in Fig. ^1 Then, the topological charge 
density is distributed in the direction perpendicular to 
the direction of polarization of the meron-pair, being con- 
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FIG. 14: (Color online) Equilibrium structures of the spin 
texture in the presence of a longitudinal magnetic field (ci) 
and a transverse magnetic field (cjr) in the region [—4 < 
x,y < +4]. The vectorial plots in the x-y plane and along 
the y=0 line are shown. The values of the parameters are 
Q. = 0.15, Co = 1000, ci = -10, C2 = 0, and (a) lu^ = 0, 
(b) LJR = 0.02, (c) t^R = 0.05, (a) ojr = 0.10. The cross 
sections of lipi^ (solid curve), \ip2\^ (dashed curve), and the 
total density px (dotted curve) along the y = line is also 
shown within the range [— 8 < a; < -|-8]. (e) The distributions 
of the topological charge density q{r) along the y = Q line 
for t^R = (sohd), 0.02 (dashed), 0.05 (dashed-dotted), 0.1 
(dotted). 



centrated on the domain-wall region. This anisotropy of 
the meron-pair gives an interesting situation when the 
condensates undergo a rapid rotation; the anisotropic in- 
teraction between the meron-pai rs g enerates a distorted 
lattice of the "vortex molecules" jl9| , which is an unique 
feature of the rotating two-component system with the 
internal coherent coupling. 



VII. CONCLUSION 

In conclusion, we have discussed the coreless vortex 
states and the corresponding spin textures in rotating 
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FIG. 15: (Color online) Equilibrium structures of the spin 
texture in the presence of a spin-spin interaction C2 and a 
transverse magnetic field (ujr). The parameter values are SI = 
0.15, CO = 1000, ci = 0, UJB. = 0.05, and (a) C2 = 20, (b) 
C2 = —20. (top) The cross sections of IV"!!^ (solid curve), |'i/'2|^ 
(dashed curve) , and the total density px (dotted curve) along 
the y — line is also shown, (middle) The corresponding spin 
texture in the region [— 4 < x, j/ < +4]. The vectorial plots in 
the x-y plane and along the y=0 line are shown, (bottom) The 
spatial distributions of the topological charge density q{r). 

two-component BECs with and without internal coher- 
ent coupling. The axisymmteric and nonaxisymmetric 
structure of the spin texture have been discussed by ex- 
ploring the NLctM derived from the Gross-Pitaevskii en- 
ergy functional. The variational function of the spin pro- 
file, which is based on the exact solution of the classical 
NLctM, provides a good description of the coreless vortex 



states. This variational method also improves greatly the 
usual analytical approach based on the Thomas-Fermi 
approximation @. In the case of SU(2) symmetry, these 
two spin textures are equivalent topological excitations 
and transform to each other by a global rotation of the 
pseudospin. We have discussed the effect of the SU(2)- 
symmetry breaking contributions, inequality of the three 
coupling constants ui, U2 and ui2 (ci and C2 term) and 
the Rabi frequency wr, on the spin textures. These con- 
tributions are interpreted as a longitudinal and trans- 
verse magnetic field and a spin-spin interaction in the 
NLcrM. For an axisymmetric skyrmion (a nonaxisymmet- 
ric meron-pair), the ci and C2 term (wr term) affect the 
optimized size of the topological excitation, but do not 
break the radial symmetry of the problem. Inclusion of 
all those terms makes a topological excitation character- 
ized by an anisotropic distribution of the vorticity and 
the topological charge density. 

One of the open problems is to study vortex states 
in a ferromagnetic regime. Rich spin textures are ex- 
pected with the combination of the ferromagnetic fea- 
ture leading to spin-domain formation and a rotational 
effect, e.g., a rotating domain-wall cross ("propellers") 
and "serpentine" vortex sheets [2^. The problem 
of their stability remains to be investigated. Further- 
more, it is necessary to study in more detail lattices of 
spin textures which appear in the rapidly rotating two- 
component BECs, where we can expect a richer phase 
diagram and structural phase transitions of the vortex 
phases 16, 43]. It is also of interest to extend the pseu- 
dospin model to two-component fermion systems in view 
of the rapid progress in ultracold fermions after the ob- 
servation of condensation of fermionic pairs |44| : rapidly 
rotating two-component fermions 45] may exhibit phe- 
nomena analogous to what has been discussed in double- 
layer quantum Hall systems. 
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